Scaling of running time of quantum adiabatic algorithm for propositional satisfiability 
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We numerically study quantum adiabatic algorithm for the propositional satisfiability. A new class 
of previously unknown hard instances is identified among random problems. We numerically find 
that the running time for such instances grows exponentially with their size. Worst case complexity 
of quantum adiabatic algorithm therefore seems to be exponential. 
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I. INTRODUCTION 

Computers play a vital role in modern society. Since 
its inception in the middle of the previous century the 
power of digital computers has been growing exponen- 
tially with time. How long can this growth be sustained? 
Great interest in recent years in quantum computing is 
partially fuelled by the discovery that quantum comput- 
ers could perform certain tasks faster than any classi- 
cal computer. An example is the famous Shor factor- 
ing algorithm which is polynomial, whereas the best 
known classical algorithm is super-polynomial. Accord- 
ing to computational complexity, problems can be di- 
vided into two large groups. Those for which the time 
to find a solution grows polynomially with the size of 
the problem belong to the so-called P class (polynomial) 
and those that require polynomial time to verify the so- 
lution belong to NP (non-deterministic polynomial). Es- 
pecially important subset of NP problems is called NPC 
(NP-complete). They have the property that any NP 
problem can be transformed to NPC problem in a poly- 
nomial time. Therefore, finding a polynomial algorithm 
for a single NPC problem would immediately provide a 
polynomial algorithm for all NP problems. Currently all 
known algorithms need exponential time to solve NPC 
problems. In a vague way it can be said that NPC are 
the hardest of NP problems. 

Recently a novel way of doing quantum computation 
via adiabatic evolution has been suggested 0, Q . The 
idea of using adiabatic evolution to do quantum com- 
putation is very simple and elegant. One starts with 
the system in the ground state of the initial Hamilto- 
nian H(0). Then the Hamiltonian is adiabatically slowly 
changed from H(0) to the final H(l), whose ground state 
encodes the solution to the problem we want to solve. 
The adiabatic theorem then ensures that if the chang- 
ing of the Hamiltonian is sufficiently slow, we end up 
in the ground state of -ff(l) at the end, thereby obtain- 
ing the solution to our problem. Numerical simulation 
of the adiabatic algorithm for a NPC problem called 
exact cover 3] indicated that quantum adiabatic algo- 
rithm might need running time that grows only quadrat- 
ically with the size of the problem. Subsequently, there 
have been many studies of quantum adiabatic algorithms, 



mostly numerical |4 | 5t l a , Ml but also some rigorous re- 
sults are known |Rl9t ll0l|llf . It has also been shown re- 
cently that adiabatic computation is polynomially equiv- 
alent to standard quantum computation |12| . Adiabatic 
algorithm has also been experimentally realized on a 
NMR quantum computer [l3| . While numerics for num- 
ber partitioning problem showed exponential scaling £| 
the scaling of running time in other NPC problems is 
still unclear. For instance, for a paradigmatic example 
of a NPC problem called 3-SAT (3-satisfiability), poly- 
nomial scaling ~ n 3 of median cost has been found |fj. 
Needless to say, the implications of having a polynomial 
quantum adiabatic algorithm for NPC problem would be 
enormous. But we have to keep in mind that the compu- 
tational complexity is defined in terms of the worst case 
performance. 

Although the studies so far focused on presumably 
hard instances of 3-SAT problems we will show that there 
exists a class of even harder 3-SAT instances not known 
before. We will present a clear numerical evidence for 
an exponential scaling of running time of quantum adia- 
batic algorithm for these 3-SAT instances. This finding 
could also be relevant for classical algorithms develop- 
ment, where hard instances are used in algorithm design 
and testing. 

The outline of the paper is as follows. In Section HT1 we 
introduce a problem studied, namely a random 3-SAT. 
In Section IIHI quantum adiabatic algorithm is defined 
and the degeneracies of the initial and final Hamiltonian 
are explored. In Section II VI we study the probability 
to successfully obtain correct result. Then in Section Ivl 
the energy spectrum during adiabatic evolution is studied 
and finally in Scction lvTl thc scaling of the energy gap and 
running time is presented. 



II. RANDOM 3-SAT 

In the present paper we will consider 3-SAT problem. 
It is a paradigmatic example of a NPC problem. 3-SAT 
formula is a logical statement involving n boolean vari- 
ables bi. It consists of m clauses Cj in conjunction (logical 
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AND = A) 

d AC 2 A---AC m , (1) 

and each clause Cj is a disjunction (logical OR = V) of 3 
literals, where a literal is a variable 6.; or its negation ^6j 
(logical NOT = -i). To illustrate, an instance of a 3-SAT 
formula with n — 4 variables and m = 2 clauses is 

(b 2 V -.63 V 64) A (61 V b 2 V -.63). (2) 

The problem is to decide whether a given 3-SAT formula 
is satisfiable, i.e. whether there exists a prescription of 
variables bi such that 3-SAT formula is true. Such pre- 
scription is called a solution, the number of which will be 
denoted by r. The formula given as an example © has 
many solutions, one being for instance (64636261) = 1101, 
where 1 and denote true and false, respectively. 

Formula that is a conjunction of disjunctions is said 
to be in a conjunctive normal form (CNF). Any logi- 
cal statement consisting of A, V and -1 operators can be 
rewritten in a CNF form. While 2-SAT problem (CNF 
formula having at most 2 literals in each clause) belongs 
to P, 3-SAT is in NPC. Besides being paradigmatic exam- 
ple of a NPC problem it has wide range of applicability in 
e.g. scheduling problems, hardware verification etc... It 
is also directly related to deductive reasoning important 
in artificial intelligence. If we are given a set E of facts 
(statements) Cj, E = UCj, a new statement C ne w can be 
deduced iff a union E U {^C ncw } is not satisfiable (we 
arrive at contradiction assuming ^C ncw ), i-e. we have to 
solve a SAT problem. 

As the computational complexity is defined in terms 
of worst case performance, hard instances of 3-SAT have 
been especially intensely studied. Frequently they are 
obtained by randomly drawing clauses, so-called random 
3-SAT. A random 3-SAT consists of m different random 
clauses, where each clause is obtained by picking 3 differ- 
ent variables bi and negating each with probability 1/2. 
This is the procedure we used to obtain random instances 
of 3-SAT. As we were mainly interested in formulae with 
only one solution, we solved each instance and rejected 
those not having exactly one solution. In the literature 
on the other hand they usually study random 3-SAT for- 
mulae with an arbitrary number of solutions. In such 
case a phase-transition is found Q with the hardest in- 
stances occurring around the phase-transition point of 
m/n ~ 4.2. As we will see, random 3-SAT instances 
with exactly one solution and small m/n will turn out to 
be harder than those at the phase-transition (and having 
arbitrary number of solutions). Correct choice of hard 
instances of 3-SAT is therefore absolutely essential in or- 
der to see clear exponential scaling of running time. For 
more information about a phase-transition in random 3- 
SAT see also papers in the 0] and also [Til Il7|. For 
classical SAT problem solving see e.g. collection Jlq. 



III. HAMILTONIAN FOR ADIABATIC 
ALGORITHM 

The construction of Hamiltonian for the 3-SAT prob- 
lem is straightforward. For each variable we have to have 
available two states, i.e. one qubit, giving the dimension 
of the total Hilbert space N = 2". The initial H(0) is 
problem independent and we choose it to be a sum of 
one-qubit Hamiltonians Hi acting on i-th qubit, 

Hi = (\ ~A . (3) 

»=1 ^ ' 

The ground state |i?o(0)) of the initial Hamiltonian has 
energy Eq = and is a uniform superposition of all com- 
putational states, 

1^0(0)) = \ b n---bi), (4) 

^ 6i,...,6„=0,l 

where a label of state |6 n 6„_i ... 61) = |b) denotes a bi- 
nary expansion of the state, i.e. the value of each qubit. 
In fact the whole energy spectrum of the initial Hamilto- 
nian is easily calculated. It consists of integer energies, 
Ei(0) — i, with the degeneracy of level Ei(0) being equal 
to (™) . The final Hamiltonian H(l) is problem depen- 
dent. We use a diagonal Hamiltonian in the computa- 
tional basis, with the energy of state |b) equal to the 
number of clauses it violates. As all states violating a 
given clause have the same values of 3 variables occur- 
ring in that clause, each clause can be represented by a 
single 3-qubit gate. For our example J2} we would have 
two terms, 

H(l) = |01G)(010| 43 2 8 1+ |100)(100| 32 i O 1, (5) 

where subscripts denote on which qubits the operator 
acts. Therefore, a state |b) satisfying all clauses (i.e. 
a solution) would have energy 0, state violating a sin- 
gle clause has energy 1 and so on. Energy spectrum 
of H(l) is therefore composed of integer values between 
-Eo(l) = for the ground state and E m (l) = m for a state 
that would violate all clauses. Of course, the degeneracy 
r of the state Eq(1) of H(l), i.e. the number of solu- 
tions, depends on the particular instance in question. In 
the thermodynamic limit, n — * 00, 3-SAT formulae with 
m/n below a "phase transition" point at m/n w 4.2 have 
many solutions, whereas formulae with m/n > 4.2 have 
no solution [l4|. Whereas the spectrum of H(0) is fixed, 
the spectrum of H(l) is instance dependent. An example 
of spectrum of H (0) and H(l) for one particular 3-SAT 
instance having n = 14 variables and m = 42 clauses 
is shown if Fig. ^ From the spectrum of H(l) we see 
for instance, that there are no states violating more than 
12 clauses and the most abundant are 3720 states vio- 
lating 5 clauses (with energy E = 5). There is only one 
state with energy E = 0, thus there is only one solution. 
Throughout the work we will focus only on 3-SAT prob- 
lems having exactly one solution, r = 1. Later we will 
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argue that such problems are expected to be the hardest. 
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FIG. 1: (Color online) Degeneracy of the initial (left) and 
final Hamiltonian (right) for one 3-SAT instance with n = 14, 
m/n — 3 and exactly one solution r = 1. Lower part of 
spectrum for all intermediate times s is shown in Fig. |2] 

The interpolating Hamiltonian between H(0) and 
H(l) is chosen according to the following prescription, 



H(s) = (1- s)H(0) + sH(l), 



t 



(6) 



where s = t/T is a dimensionless time and T is the total 
running time of the adiabatic algorithm. The evolution 
of an arbitrary state is given by a time dependent 
Schrodinger equation 



(7) 



where for our choice of constant speed along the interpo- 
lating path H(s) (dt/ds) is constant and equal to T. 
We set H = 1 throughout the paper. 

The interpolating "path" between H(0) and H(l) is 
quite arbitrary as well as the initial Hamiltonian H(0), 
while the final Hamiltonian is determined by the prob- 
lem in question. Instead of having uniform speed of in- 
terpolation we could vary it according to the energy gap. 
Such refinements are not the subject of the present paper. 
Trough adiabatic quantum computation hard question of 
computational complexity is translated into (perhaps?) 



easier question of the scaling of the energy gap. At least 
there are plenty of tools available for studying energy 
gaps. 



IV. FAILURE PROBABILITY 

In studying quantum adiabatic algorithms one usually 
numerically looks at the probability of successfully find- 
ing the solution for different running times T l&S.blul- 
The adiabaticity condition, guaranteeing adiabatic evo- 
lution, is usually stated as 



\\dH/ds\ 



ds, 



(8) 



where j(s) — Ei(s) — Eq(s) is the energy gap between 
the ground state and the first excited state. Adiabatic 
condition can also be stated locally, saying that the local 
inverse speed has to be larger than 



dt ^ h \\dH/ds\ 
ds j{s) 2 



(9) 



The above two adiabatic conditions JSJ give us the 
necessary condition for the adiabatic evolution. What 
would be desirable to know is also the probability of non- 
adiabatic transitions. This would then give us a direct 
way to calculate necessary running time for the desired 
probability to stay in the ground state at the end. For 
2 and 3 level systems there is an exact expression. In 
the limit of slow evolution it is the famous Landau-Zcncr 
formula |l9l l20| , giving the probability of a transition in 
a 2 level system, where the two eigenenergies are E 0> i (s), 
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(10) 



here A is a minimum gap 7(s), T is a constant parame- 
ter connecting s — t/T and t runs from — oo to oo. For 
a discussion of transitions in multilevel systems see |2lj . 
Landau-Zener formula has been used before to de- 
scribe adiabatic algorithm under the assumption that 
random matrix statistics applies to avoided crossings. 
Two level transition probability is also used in the adia- 
batic algorithm for Grover search [9| . 

Let us first have a look at the lower part of a typical 
spectrum of H(s) to see what is the nature of transi- 
tions. To find the lowest eigenvalues we used implicitly 
restarted Lanczos method [22|, suitable for sparse eigen- 
value problems. By this we could find lowest few eigen- 
values for n upto 20, i.e. Hilbert space sizes N ~ 10 6 . 
One particular example for n — 14 and m = 42 is shown 
in Fig. [2 Only lowest 18 levels out of total N = 16384 
are shown. We can see that there is only one avoided 
crossing, the same was the case in all other cases we have 
checked. The reason to have only one avoided crossing 
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FIG. 2: (Color online) Lower part of spectrum (lowest 18 
levels out of total N = 16384) for one instance of 3-SAT with 
n = 14, m/n = 3, the same instance as in Fig. We can 
see avoided crossing at s ~ 0.7. Dashed line is fitted avoided 
crossing, i.e. line -Ei(s) — \J A 2 + (As) 2 , whose parameters 
are obtained by fitting parabola at the crossing. 



S 0.1 

Q. 



0.01 



10 



100 




10000 



FIG. 3: (Color online) Comparison of success probability 
p(ground) obtained with direct numerical simulation (pluses) 
with Pf obtained from Landau- Zener formula ill (H . In the 
inset the same data is shown in a semi-log scale for longer 
times. All is for the same 3-SAT instance shown in Fig. |5] 



out the calculation. At the end of the simulation, at 
time t = T, we obtain a final state \ijj(T)) and then cal- 
culate the overlap with the ground state, giving us the 
numerical probability p(ground) to remain in the ground 
state. In Fig. [3]we compare this numerical value with the 
Landau-Zener formula. One can see a very good agree- 
ment already in a regime of small T, where the probabil- 
ity to stay in the ground state is small. Therefore, the 
Landau-Zener formula perfectly describes the probabil- 
ity to stay in the ground state in all practically relevant 
regime (i.e. for high p(ground)). From now on we will 
focus on the scaling of tlz with the size n of a problem. 



is unclear to us. It might be connected with the fact 
that we have a finite gap at the beginning and at the 
end of the algorithm. Similar behaviour has been found 
by other researchers |6j . Now provided we have only one 
close encounter of two lowest levels, we can use degen- 
erate perturbation theory, i.e. consider only two closest 
levels. Transition probability is then given simply by the 
Landau-Zener formula pO|l . We first checked how well 
the Landau-Zener formula describes the transition prob- 
ability for a real 3-SAT case, such as is for example the 
one in Fig. EI To obtain probability for a transition from 
the ground state using Landau-Zener formula (|10fl we 
fitted parabola around the avoided crossing to determine 
two parameters, the gap A and the asymptotic slope A, 
needed in the Landau-Zener formula. For a case in Fig. [3 
for instance, we obtain A = 0.0318 and A — 2.67 (dashed 
line in Fig.|2J). Then we compared Pj (|10fl with the re- 
sult of a direct numerical simulation of a time dependent 
Schrodinger equation @ . We discretized time into small 
steps dt (typically dt ~ 0.1) and then calculated one-step 
propagator U(dt) = exp (—iH(t)dt/h) by expansion in a 
power series. The precision has been controlled through- 



V. SPECTRUM 

There are two parameters determining tlz- The 
asymptotic slope A at the avoided crossing has value of 
around 2 (ground state energy changes from to w 1 at 
the avoided crossing and back to at s = 1) and does not 
vary appreciably with n. The main dependence of tlz on 
n will be therefore given by the scaling of the minimal 
gap A. Heuristically one could argue that the gap A will 
be smaller when more excited levels are crowded into a 
region of energies between E « 1 and E«2, where the 
avoided crossing takes place. The number of such lev- 
els is connected with the degeneracies of the first excited 
levels at the beginning and at the end of the algorithm. 
The degeneracy of the first excited level of H (0) is fixed 
and equal to n, whereas the degeneracy of the first ex- 
cited level at the end varies from instance to instance. 
We therefore first wanted to identify instances that have 
the highest number of first excited states at the end s = 1 
(remember that we always look at 3-SAT problems hav- 
ing exactly one solution, r = 1). We calculated average 
degeneracies (averaged over 10000 instances of random 
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3-SAT with r = 1 solution) of the final hamiltonian H (1) 
for n — 10 and different m. The results are in Fig. 21 
We can see that the peak of the degeneracy moves to 
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FIG. 4: (Color online) Degeneracies of the spectrum of final 
hamiltonian H(l) for different ratios m/n. Degeneracy of 
the first excited state increases with decreasing m/n. All 
points are an average over 10000 random 3-SAT instances 
with exactly one solution, i.e. the degeneracy of the ground 
state with E = is r = 1. 

smaller energies with decreasing m/n. More importantly, 
the degeneracy of the first excited state (the height of 
the second bar) increases with decreasing m/n as well. 
Therefore, to have as high degeneracy of the first excited 
state as possible, we have to have small ratio m/n. Note 
that nothing particular happens at the point of "phase- 
transition" at m/n « 4 — 5. From now on we will choose 
m/n = 3, as such 3-SAT instances have high degener- 
acy of the first excited state. If this degeneracy increases 
exponentially with n we have a fair chance that the gap 
A will also decrease exponentially. We therefore checked 
the scaling of the degeneracy of the first excited state of 
H(l) with n at a constant ratio m/n = 3. Results are 
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FIG. 5: (Color online) Degeneracy of the first excited state of 
the final -ff(l) for different n and m/n. Exponential growth 
can be seen for m/n — 3, while for m/n = 5 and 4.5 (around 
the phase transition) the degeneracy is much smaller and the 
functional form of the growth can not be established. 

shown in Fig. [SJ We can see clear exponential growth 



for m/n = 3 with the fitted exponential « 8.2e 021 ™. For 
larger m/n, i.e. around the phase transition, the growth 
is much slower and an exponential dependence can not 
be firmly established. 

Random 3-SAT instances with m/n = 3 and exactly 
one solution are therefore our candidates for hard prob- 
lems. For such instances the number of states violating 
only one clause (i.e. the degeneracy of the first excited 
state of H (1)) grows exponentially with n and so the en- 
ergy gap is expected to decrease exponentially and con- 
sequently running time to increase exponentially. What 
is the difficulty of such 3-SAT instances for classical algo- 
rithms? One could expect that classical algorithms will 
also have a hard time finding the solution as it is "sur- 
rounded" by exponentially many states that violate only 
one clause. To illustrate this we tested classical GSAT al- 
gorithm with random walk extension [2j| . GSAT belongs 
to a group of incomplete algorithms for 3-SAT meaning 
that there is no a priori terminating condition. All such 
algorithms are variants of a local search. Concretely, in 
GSAT algorithm one variable is negated (flipped) at each 
step. The variable to negate is chosen so that the result- 
ing state satisfies a maximum number of clauses. Note 
that GSAT algorithm can not prove unsatisfiability and 
is therefore suitable only for solvable instances. In Fig. El 
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FIG. 6: (Color online) Running time (number of flips) for 
classical GSAT algorithm solving 3-SAT instances with ex- 
actly one solution (top points) and instances with at least 
one solution (bottom points). All is for n — 20 and an aver- 
age over 1000 random 3-SAT instances is performed. Error 
bars show standard deviation at each point. 

we show the average number of flips needed to find a solu- 
tion for n = 20 and different number of clauses m. Data 
shown is an average over 1000 instances having exactly 
one solution, r = 1, and for comparison also for problems 
with at least one solution, r > 1. We can see, that for 
problems with at least one solution we have a characteris- 
tic phase-transition dependence [l4l Il5j with the hardest 
instances at around m/n w 4.2. On the other hand, for 
r = 1 the difficulty of problems grows with decreasing 
m/n. We see that such problems are actually harder 
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than the problems at the phase transition. The same be- 
haviour is expected also in other local search algorithms. 
On the other hand there are also complete methods for 
solving satisfiability problems. The most widely used is 
the so-called DPLL 24] algorithm and its derivations. It 
searches trough a solution tree and can both prove un- 
satisfiability or find a solution in a finite (exponential) 
number of steps. Preliminary results show |25j that ran- 
dom 3-SAT instances with one solution are harder than 
instances at the phase-transition also for such complete 
algorithms. 
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FIG. 7: (Color online) Frequency of 3-SAT formulas, i.e. the 
inverse probability to get such a formula, with exactly r = 1 
solution among random instances at m/n — 3. Each point is 
an average over 100 instances. Line is an exponential fit. 

Why have been such instances with one solution and 
small m/n overlooked so far in a vast literature on phase- 
transition in random 3-SAT 0EE E3 ? The answer 
is very simple. They are exponentially rare among ran- 
dom 3-SAT instances (for small m/n most have many 
solutions) and so they do not show up in the average be- 
haviour that was usually studied. Nevertheless, as the 
computational complexity is defined in terms of a worst 
case performance, such problems are important. In Fig.[7| 
we show how frequently one gets a 3-SAT problem with 
exactly one solution among randomly drawn 3-SAT prob- 
lems at small m/n — 3. The best fitting line in the fig- 
ure is exponential w 2.3e 4 ". The fact that 3-SAT with 
r = 1 and small m/n are exponentially rare makes their 
generation very time consuming as we have to solve very 
many instances before we arrive at the "right one" with 
exactly one solution. 



VI. SCALING OF RUNNING TIME 

Finally, after identifying 3-SAT instances that are ex- 
pected to be hard for quantum adiabatic algorithm due to 
exponentially many states residing just above the ground 
state and after seeing that such instances are hard also for 
classical algorithms, we turn to the numerical calculation 



of the scaling of running time of quantum adiabatic algo- 
rithm for 3-SAT problems with r = 1 and m/n — 3. For 
each n we generate 100 such instances, find the position 
of the minimum of the energy gap and from the curvature 
at the gap determine the necessary running time accord- 
ing to the Landau-Zener formula (|10|) . In Fig. [S] we show 
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FIG. 8: (Color online) Dependence of the energy gap A on n 
for 3-SAT with m/n = 3 and r = 1 solution. Top points are 
for the average gap and bottom for minimal gap (average over 
100 instances). Two lines are exponential fits. For n = 16 and 
n = 18 we show the results of two independent runs to give 
the impression about fluctuations. 

the dependence of the minimal gap A on the size n. Ex- 
ponential fit gives dependence w 0.9e~ 02n for the average 
A and ~ 2e~ 52 ™ for a minimal A (out of 100 instances). 
Clear exponential decrease can be seen over several or- 
ders of magnitude. Similar exponential behaviour can be 
seen also in the dependence of the running time tlz in 
Fig. Exponential fits give scaling « 1.4e 0,9n for a max- 
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FIG. 9: (Color online) Dependence of running time (according 
to Landau-Zener formula) on n for 3-SAT with m/n = 3 and 
r — 1 solution (same data as for Fig. |HJ • Top points are for 
a maximal gap, middle for the average and bottom for the 
median tlz (average over 100 instances) . Exponential growth 
over 7 orders of magnitude can be seen. 
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imal tlz, ~ 0.05e°- 9ra for the average tlz and « 0.2e ' 6 ™ 
for a median time. We have therefore numerically estab- 
lished exponential growth of running time of the adiabatic 
quantum algorithm for 3-SAT with the problem size. 

We should mention that if one looks at the scaling of 
e.g. energy gap for problems with one solution at larger 
m/n, say around the phase-transition, the functional de- 
pendence can not be firmly established. It looks like a 
power-law ~ n -0 8 , but the exponential dependence with 
a very small exponent (e.g. ~ e -o.05n-j canno ^ \y e ex _ 
eluded. This is in agreement with a similar behaviour 
of degeneracies in Fig. \5\ It might be that the asymp- 
totic regime of exponential dependence is not yet reached 
for small n < 20 amenable to numerical study. This is 
probably the reason why in previous studies no clear ex- 
ponential time dependence could be identified. 

VII. CONCLUSION 

We have numerically studied quantum adiabatic algo- 
rithm for 3-SAT problem. First, we identified a class of 
difficult random 3-SAT instances not known previously. 
These are instances with exactly one solution and small 
number of clauses (e.g. m/n — 3). Such 3-SAT instances 
are exponentially rare, which is the reason they have not 
been observed so far. Nevertheless, as the computational 
complexity is concerned with the worst case performance, 
they determine the complexity of the algorithm. They 



have exponentially many assignments of variables that 
violate only one clause and thereby exponentially many 
states just above the ground state of the final Hamilto- 
nian. Therefore, the energy gap for such instances de- 
creases exponentially with the size of the problem and as 
a consequence, running time grows exponentially. This 
provides a firm numerical evidence that the usual quan- 
tum adiabatic algorithm for 3-SAT has exponential com- 
plexity. 

In addition, such class of 3-SAT instances is expected 
to be difficult also for classical algorithms. In local search 
algorithms exponentially many "fake" solutions, violat- 
ing only one clause, will effectively "shadow out" the real 
solution. In complete methods, like DPLL, a wrong as- 
signment of a variable early in the search tree will cause 
large backtrackings. Becouse problems with small m/n 
are under constrained, such wrong assignment will be 
very frequent which will make the search tree very large. 
As the performance of 3-SAT solving algorithms is im- 
portant in many areas, further study of the behaviour of 
classical algorithms on this new class of 3-SAT problems 
in necessary. 
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